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High-sensitivity imaging with multi-mode twin beams 
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Twin entangled beams produced by single-pass parametric down-conversion (PDC) offer the op- 
portunity to detect weak amount of absorption with an improved sensitivity with respect to standard 
techniques which make use of classical light sources. We propose a differential measurement scheme 
which exploits the spatial quantum correlation of type II PDC to image a weak amplitude object 
with a sensitivity beyond the standard quantum limit imposed by shot-noise. 



I. INTRODUCTION 
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2 ', The field of Quantum Imaging aims to exploit the quantum nature of light and the natural parallelism of optical 
' signals to devise novel techniques for optical imaging and for parallel information processing at the quantum level. 
_^ , To cite a few examples, multi-mode quantum correlations in macroscopic tvi^in beams have been used to measure 

■ small laser beam displacements beyond the Rayleigh limit [3| and for the noiseless amplification of optical images [2] ■ 
Parallely, a number of new applications that use multi-mode twin beams in the low gain regime have been proposed, 

I— I' such as e.g. quantum lithography |3jj|, entangled photon microscopy ^Jj ^ind dispersion canceled quantum optical 
\ coherence tomography Q . The state of the art in this field can be appreciated by reading i . An old review is given 
II ' by 0], while a new review is in press Q. 

In this paper we focus our attention on a specific application, namely the possibility of improving the sensibility in 
the measurements of very weak images, i.e. the intensity distribution transmitted by objects with a small absorption 

■ coefficient. The detection of a weak amount of absorption (but not of its spatial distribution) with a sensitivity beyond 
the standard quantum limit (SQL) was demonstrated in the past by using single-mode twin beams produced by cw 
optical parametric oscillators (OPOs) [§, [13, [lHj [l^l ■ these experiments the sub-shot- noise intensity correlation 

I [ of the twin beams was exploited in order to perform differential intensity measurements of a very slight amount of 
^ . absorption. For example in this technique was used in order to measure a slight amount of absorption from an 
' electro-optical amplitude modulator, achieving 7dB of noise reduction with respect to the SQL of a classical differential 
lO ■ scheme. In pd| a spectroscopic measurement of a two-photon transition was implemented with a noise reduction of 
O ■ 1.9 dB. 

] However, single-mode twin beams cannot be exploited to retrieve information on the spatial distribution of the 
• transmitted field since the correlation vanishes as soon as one detects small portions of the two beams instead of 
T-H ■ the whole beams. Being interested in measuring an image, we are forced to consider a multi-mode source which is 
C . able to display quantum correlation also in the spatial domain. Recently, our group demonstrated the existence of 

■ such kind of correlations in the high gain regime of single-pass parametric down conversion (PDC) BlilQIIl. In 
J> , the experiment [3, the sub-shot noise correlation between symmetrical points of the PDC far field was observed 

and was interpreted as a manifestation at the macroscopic level (i.e. in the large photon number regime) of the 
^ . , transverse momentum conservation of photons. The aim of this paper is to show that this spatial twin beam effect 

■ existing over the several phase conjugate signal and idler mode pairs offers the opportunity to retrieve the full 2D 
spatial distribution of a weak object in parallel, with an improved signal-to-noisc ratio with respect to standard 
techniques. This opportunity is interesting, for example in the case of biological samples or whenever there is the 
need for illuminating the object at low intensities. 

The paper is organized in the following manner. In Sec|TT]we briefly review the principle of operation of differential 
measurements used to detect faint amounts of absorption, showing with a simple two-mode model that the use of a 
source displaying quantum correlations can provide a higher sensitivity than a classical source. In Sec lIIII we propose 
an imaging scheme that can be used to detect the spatial distribution of amplitude objects with a PDC source. The 
model equations that describe the multi-mode PDC process are illustrated in Sec lIVI In the last sections of the paper 
(Sec|V]to Sec lVIII) we investigate on the robustness of the spatial quantum correlation necessary to implement the high 
sensitivity imaging technique, considering experimental imperfections that are difficult to avoid in a real experiment. 
In particular in Sec. IVIII we present the results obtained through a numerical stochastic model that simulates the 
imaging experiment with realistic parameters. 
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II. DESCRIPTION OF THE DIFFERENTIAL DETECTION SCHEME 



A common procedme to detect a weak object makes use of a differential measurement. Tliis technique is illustrated 
schematically in Fig. [1] In its classical version [Fig. [TJa)] a laser beam (not necessarily shot noise limited) is separated 
with a 50/50 beam splitter into two classical "twin" beams, a "test" beam which illuminates the object under 
examination, and a "reference" beam which does not interact with the object. In the ideal case in which the two 
beams are perfectly balanced, they can be considered as two classical copies one of each other, meaning that their 
intensity fluctuations display the strongest level of correlation allowed by a classical source (corresponding to the 
shot noise level of the incident beams). The subtraction of the detected intensities allows to eliminate the classical 
excess noise contained in the source field and to retrieve information about the object with a far better signal-to-noise 
ratio (SNR) than through direct illumination. However, when a classical source is used, this technique is limited 
intrinsically by the level of shot noise. 

We shall now illustrate how, under appropriate conditions, quantum correlations allow to achieve a better SNR of 
the object than classical correlations. In the quantum scheme [see Fig. [ijb)], the reference and the test beams are 
replaced by the signal and idler fields generated by a PDC source, for example the single-mode twin beams generated 
by an OPO above threshold. 
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FIG. 1; Scheme for the detection of a weak object through a differential measurement. In (a) the test and the reference beams 
are obtained by splitting a laser beam with a symmetric beam splitter; in (b) the signal and idler beams of a PDC source are 
used. 



Let us indicate with ai and a2 the field operators of the two beams at the plane before the object, and with 
Nj — ajcj, j — 1,2, the corresponding photon fluxes. The fields Ui and 02 obey the usual canonical commutation 
relations 

[a,, a]] = , = 1,2). (1) 

The degree of correlation of the two beams is determined by the ratio between the variance of iV_ = N2 — Ni and the 
corresponding level of shot-noise, i.e. 

(iVi) + (iV2) ' ^ ' 

It is equal to unity in the case of the classical copies obtained with a well balanced 50/50 beam splitter, while it can 
go well below unity when the PDC source is employed (it vanishes in the ideal limit of perfect quantum intensity 
correlations). We assume that the two beams are perfectly balanced and have symmetric statistical properties, so 
that in particular 

(TVi) = {N2) , (SN!) ^ (5Ni) . (3) 

Indicating with tobj the transmission coefficient of the object, the field in the test arm undergoes the unitary trans- 
formation 



a'l =tobj ai+i^l-tl^^v , (4) 
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where v denotes here a vacuum field operator. The measured quantity is the difference of the photon fluxes collected on 
the two detectors after the test beam has passed through the object, associated with the operator N'_ = a\a2 — a\ 'a'^. 
Recalling the condition (A^i) = (-/V2), one immediately sees that the mean value of N'_ is proportional to the object 
absorption coefficient a = 1 — |iobjp: 

{N'_) - (4a2> - |tobjP(aI«i> = (^i>« (5) 

The noise of the measurement is determined by the variance of N'_ . By using relations ((!]) and ^ together with the 
symmetry conditions ([3]), we can express {5N']:) in terms of the unprimed quantities (iVi), {SNf) and {6N'^), which 
depend on the photon statistics and correlation of the two beams in absence of the object: 

{SN'_^) = a^{6Nl) + (1 - a){6Nl) + a{l - a){Ni) . (6) 

Using ([6]) together with definition ([2]) we obtain the SNR of the measurement 



SNR.. <£^ = ^^^^= (7) 

'{SN'^) \/a'^En + 2a{l- a) + a 



where we have introduced the quantity 



- — m — 



which characterizes the signal (idler) noise in excess with respect to that of a coherent beam with the same intensity 
(it is the so called Mandel Q-factor, used when dealing with sub-poissonian statistics for which Q takes negative 
values). 

In the classical scheme of Fig. [IJa), provided the two beams are perfectly balanced and satisfy conditions ([3]), it 
can be verified that a — 1 (i.e. {SN"^) — 2{Ni)) whatever the amount of excess noise. Thus we find 



SNR.,... = ^ . (9) 

Va^E„ + 2 -a ^ ' 

The standard quantum limit is obtained from this expression by considering the case of a coherent beam, for which 

En = and 



SNR,,, ^ . (10) 

Noting that for a weak object the condition En « is easy to achieve, the SNR obtained with the differential 
measurement ^ usually approaches the SQL value PH)) even if the classical source is not shot-noise limited but 
presents some excess noise. 

It is worth noting that the SQL could also be defined by considering the SNR of a direct measurement scheme, which 
makes use of a single coherent beam with the same mean photon number (iVi): this would give SNRgg, = a^J (iVi), an 
improvement by a factor \/2 with respect to the SQL of the differential measurement (jlOp. However in the following 
we will use definition (fTO|) since the detection of weak amount absorption is most commonly obtained through the 
differential measurement technique that allows to approach the SQL much more easily than direct measurements. 

Let us now investigate the conditions under which the SQL can be beaten by using photon number correlated 
beams. To this end we introduce the ratio 



SNR^^ 2-a 
SNRsQL Y a2S„ + 2(t(1 - a) + a ^ ' 

which measures the improvement of the SNR achieved by using the PDC source with respect to the SQL. In the limit 
of a weak object we have 

R^-^ ior a «1, a^En^2a (12) 
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Thus, provided the excess noise is not too large, the SNR is improved by a factor ^ I/a/o' with respect to the SQL. 
More generally, it is easily seen that R > 1 provided that 

o- < (TMAX = 1 - r (13) 

Z[l — a) 

fJMAX represents the maximum value of a for which we have an improvement in the SNR with respect to the SQL. 
Clearly we have always omax < 1 (since i?„ > 1 and < a < 1). This is an expected result: the twin beams of 
the PDC source need to be correlated below the shot noise level (i.e. cr < 1) in order to beat the classical source 
configuration. It can be seen that condition (jl3p is easily verified when a — s- since Umax is close to unity in this 
limit. However the goal is to obtain a substantial improvement of the SNR and Eq. shows that the condition on 
a is rather severe. 



III. THE HIGH-SENSITIVITY IMAGING SCHEME 



The previous treatment did not consider the spatial aspects of the detection process. Since we are interested in 
developing an imaging system, parallel multi-mode operation is required. We focus here on the nonclassical source 
case, being interested in determining the conditions where the SQL can be beaten. The process of spontaneous 
parametric down-conversion is naturally broadband in the spatial as well as in the temporal domains and is therefore 
well suited for this task. In particular, we can exploit the sub-shot noisephoton-number correlation arising between 
pairs of symmetrical regions of the far field, theoretically predicted in [ij, [l^ and experimentally demonstrated in 
[l4i [Tsj . Close to the degenerate frequency, photon pairs are indeed emitted almost symmetrically with respect to the 
pump axis because of the conservation of the transverse momentum. The finite waist Wp of the pump beam introduces 
a spread in the relative angular directions of the twin photons on the order of the pump angular bandwidth cx X/wp, 
A denoting the signal/idler wavelength. Assuming the far field is observed in the focal plane of a lens with focal length 
/ (see Fig. [2]) this angular spread corresponds to the coherence length of the field in this plane Xcoh ~ Xf /wp. In 
it was shown that the fluctuations in the photon number difference detected in two symmetrical regions go well below 
shot-noise (ideally to zero when losses are negligible), provided the detection areas are large compared to Xcoh- 




FIG. 2: High sensitivity imaging scheme that exploits sub-shot noise correlations between symmetrical angular directions of the 
PDC field. The far field is observed in the focal plane of a lens with an /-/ imaging system. The object absorption coefiicient 
a{x) is retrieved by subtracting the intensities measured from symmetrical pixels of the two arrays of detectors. 

In accordance with this result, we consider the imaging scheme illustrated in Fig. [2] The pump field which illuminates 
the x^^'' nonlinear crystal is a coherent pump pulse with a large beam waist Wp. The object, characterized by an 
absorption coefficient a{x) depending on the transverse plane position x, is placed in the far field of the source. The 
lens shown in the figure is used to perform the Fourier transformation of the fields from the source into its second focal 
plane, which is taken as the detection plane. Both the signal (test) and the idler (reference) photons are collected by 
two arrays of detectors, typically the pixels of a high quantum efficiency CCD camera [13] • The detectors are either 
placed immediately after the object plane as in the figure, or more realistically, a telescopic system is used to image 
the object plane into the detection plane. The far field coherence length Xcoh ~ Xf /wp is assumed to be on the same 
order of magnitude of the pixel detector size or smaller, in order to fulfill the sub-shot noise correlation condition and 
compatibly with the requirement of large photon number operation. In this way, the ffuctuations of the reference and 
the test beams display strong spatial {point to point) correlation for each symmetrical pair of detectors of the CCD 
array. Sub-shot noise correlation for each of these pairs of detection channels is the goal to achieve in order to beat 
the SQL in a differential measurement. 
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IV. THE NUMERICAL MODEL 



With the aim of simulating the high-sensitivity imaging experiment proposed in the previous section, we developed 
a fully 3D numerical model which takes into account realistic features of the PDC process, such as the phase-matching 
conditions inside the crystal as well as the finite size of the pump beam both in the temporal and the spatial domains. 
The model is based on a set of equations describing the propagation of the slowly varying envelope operators of the 
signal and idler fields, which we denote by ai{z, x, t) and 02(2, x, t), and satisfy boson-like commutation rules at equal 
z, i.e. [ai(z, af, t), aj(z, X, i)] = Sij6{x — x')6{t — t'), = 1,2). As we are considering the parametric regime, it is 
assumed that the pump field Ap(x,t) remains almost undepleted and can be treated as a known classical field. For 
definiteness, we shall assume it has a Gaussian profile with a plane-wave front at 2: = 0, with a beam waist Wp and 
duration Tpi 

Ao{z = 0, X, t) = (27r)-3/2^pe-(^'+«')Ke-*'/^p . (14) 
The propagation equation for the signal field oi can then be written in Fourier space as (author?) [TgI. [T§|: 
d 

— ai(z,g, ri) = I [fci2(g, r2) - fci]ai(z,9, r2) (15) 



dz 



I 9_ -iAaz 
It. 



the equation for the idler wave 02 can be obtained by exchanging the subscript 1^2. The operators 0^(2,(7', 51) 
denotes the Fourier transformation of the fields aj{z,x,t): 

a,(z, q,n)^J^J ^a,{z, x!, t)e-?-^+^"' (j ^1,2). (16) 

where f2 denotes the temporal frequency relative to the carrier frequency ujj of field j, q the transverse wave vector 
component; kjz(q,i^) = yjkj{q, 51) — q^ is the corresponding longitudinal component along the propagation direction 
z. The wave-number function kj{q^ il) determines the linear propagation properties of field j; in particular its linear 
and quadratic dependence on f2 is responsible of temporal walk-off and second-order temporal dispersion, while its 
linear and quadratic dependence on g leads to spatial walk-off and diffraction. The parametric gain 17 is a dimensionless 
constant proportional to the strength of the nonlinear interaction, to the crystal length Ic, and to the peak value of 
the pump field Ap] its value determines the number of photons that are generated in the down-conversion process in 
mode pairs that are well-phase-matched. Aq = fci + fc2 — fcp indicates the collinear phase-mismatch parameter of the 
three carrier waves. The normalized function ap((f, il) = 1/ {tt^^'^ 5q^j5ujp)e~'' l^ipe~^ l^^p is the Fourier transform of 
Aq(x^€) I Ap^ with bqp = ^jwp and bujp = 2/rp denoting the pump spatial and temporal bandwidths respectively. 

In the plane- wave pump (PWP) limit where bqp and bujp are much smaller than the spatial and temporal bandwidths 
of phase-matching, Q!p((f, 51) 5((j)5(51) and the model equations can be solved analytically. The solution of Eqs. (fT51) 
can be written as an input-output relation in the form of a two-mode squeezing transformation: 

ar* (g, f^) = Ui (g, f^)ar* (9, + ^1 (9, ^^)4" (-9, , 

aT\q, 51) - [/i(q, n)ar\q, n) + V,{q, 11)4™ (-g, > (17) 

The explicit expressions of the gain coefficient Uj can be found in [l^, [l^. Relation pT|) links only phase-conjugate 
modes {q, 51)i and {—q, —51)2 of the signal and idler fields. This feature reflects the perfect correlation of the transverse 
momenta and energies of the emitted twin photons which occurs in the PWP limit. In particular, the signal-idler 
spatial correlation function is found to display a delta-like peak in correspondence to opposite positions in the far 
field plane, as a consequence of the transverse momentum conservation rule [Tsl. flGj). 

As already mentioned in the previous section, the finite pump bandwidth Sqp in the spatial domain introduces an 
indeterminacy in the correlation between the angular directions of the emitted twin photons. This indeterminacy 
gives rise to a spread of the signal-idler correlation function which acquires a finite width on the order of the far field 
coherence length Xcoh ~ /wp- Thus, as a consequence of the finite transverse size of the pump beam, perfect spatial 
correlations between two symmetrical ideal detectors can be achieved only if the linear size of the detectors is large 
compared to Xcoh] in this ideal limit we would have a — 0. Conversely, if this condition is not fulfilled the fluctuations 
of 7V_ just lie below or close to the shot-noise level, i.e. we have < <t < 1. When the flnite quantum efficiency of the 
detectors 77 < 1 is taken into account, the highest level of correlation that can be reached is cr = 1 — 77 (see e.g. jT3|). 

In the simulations we shall consider an optical setup close to that of the experiment performed in Como to measure 
sub-shot-noise far field correlations in the high gain regime of PDC (author?) [T^i [3 '■ the nonlinear crystal is a 4 mm 
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long beta-barium borate (BBO) crystal cut for type II phase-matching and PDC is observed around the degenerate 
wavelength A — 704 nm. The far field is measured in the focal plane of a lens of focal length / = 5cm as shown in 
Fig. [21 The pixel detectors of the CCD are characterized by a 20x20 fiui^ square area and a high quantum efficiency 
77 approaching 0.9 at the signal wavelength. The effective quantum efficiency which takes into account the losses of 
all optical elements was evaluated to be around 0.75 [T3 |. 

All the numerical results presented in the following sections have been obtained by considering a stochastic model 
based on the Wigner representation equivalent to Eqs. ([T5)) . All the expectation values of the field operators moments 
necessary to evaluate the observables of interest were evaluated by integrating numerically the classical-looking prop- 
agation equations of the model (formally identical to Eqs. (fTS)) ) and by performing the necessary stochastic averages 
on the output fields. As the Wigner representation provides only symmetri cally ordered operator moments, we had 
to apply appropriate corrections to obtain the desired ordering (see e.g. [H, [3 more details). 



V. BEHAVIOUR OF THE SIGNAL-IDLER FAR FIELD CORRELATION FUNCTION 



In this section we briefly discuss some aspects of the signal-idler intensity correlation that are relevant for the 
proposed imaging scheme. The spread of this function indeed determines the minimum size of the detection area for 
which it is possible to observe sub-shot noise correlation and sets therefore a lower limit for the imaging resolution. 
To be more precise, let us consider the normalized correlation function for the signal and idler intensity fluctuations, 
defined as 

1 i2(a;i, a;2 ) = — , (18) 

where SNj{xj) = Nj{xj) — {Nj{xj)), {j — 1,2) denotes the photon number fluctuation operator in position Xj of the 
detection plane for field j. Notice that with this normalization |ri2| < 1 and = 1 represents the maximum possible 
amount of correlation. As discussed in the previous section, the correlation function displays a peak for X2 = —xi as 
a consequence of the q —q momentum correlation of the emitted twin photons. 
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FIG. 3: Cross-section of the normalized signal-idler correlation function ri2(a'i,— a;i -|- x) as evaluated from the stochastic 
numerical model for increasing values of g. The pump beam waist is Wp — 1500//m, the lens focal length / = 5cm. The spatial 
step size of the numerical grid is bfim. 



At low gains, i.e. for 5 << 1, we can use a first order perturbative expansion of the solution of Eqs. (|15p in power of 
g in order to evaluate ri2(a:i, 0^2 )• It is found [20j | that it is proportional to the square modulus of the pump Fourier 
transform ap{q, fl) evaluated in (aJi -I- a!2)A//27r. According to this result, and considering the Gaussian pump proffie 
of Eq. (HH), the FWHM of the correlation function is then given by 



_ V2ffi2 A/ 

xpwHM — , Uyj 

TT Wr, 
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a value corresponding to the typical far field coherence length Xcoh- With the parameter values / = 5cto, A = 0.704^m 
and Wp — 1500/ini, we obtain xfwhm = 8.8/ini. 

For values of g on the order of unity or higher, no analytical expression of T12 is known and we have to resort to 
the numerical stochastic model described in the previous section. Fig. [3] shows the result of these simulations: the 
cross-section of Ti2{xi^—xi + x) discretized on the numerical grid is plotted as a function oi x = {x,y) along the 
X-axis (corresponding to the walk-off direction) for increasing values of the gain parameter g. It has been evaluated 
by performing ensemble averages as well as spatial averages, the latter being allowed by the translational invariance 
of the signal and idler fields in the considered far field region. In order to have a significant resolution, here and in 
the results that follows we consider a numerical grid with a spatial step size of 5/im in the far field plane (according 
to the mapping q <-> (A//27r)q* holding between the spatial frequency plane and the focal plane of the /-/ lens system 
illustrated in Fig[21). As can be seen from the figure legend, the smallest value g — 1 yields a FWHM of ~ 11/im, close 
to the low gain limit result (jl9p . However, by increasing the parametric gain we observe that the correlation function 
becomes broader and broader: its FWHM is more than doubled with respect to the low gain value as g is raised to 
4.5. This stronger spread occurring at high gains can be explained with the following intuitive arguments. Inside the 
crystal, the cascading effect which causes the exponential grow of the number of generated photon pairs is enhanced 
in the regions where the pump field takes its highest values. Thus, in a regime of very high gain, most of the photon 
pairs are produced where the pump field is close to its peak value, i.e close to the center of the beam. As a result the 
effective region of amplification inside the crystal becomes narrower than the pump beam profile (see also [Tsl. [2]|). 
and this effect produces a broadening of the signal-idler intensity correlation function in the far field. 

Finally we note that the peak value Fy^"^ of the discretized correlation function is related to the the correlation factor 
between symmetrical pixels, a, and the excess noise on those pixels, E^, through the relation F^"^ = 1 — a/{En + 1); 
this explains the behaviour of F^°^ which approaches unity as g (and hence En) is increased (see Fig 13]), although 
the degree of correlation (measured by a) usually deteriorates at high gains, as we shall show in the next section. 



VI. FRAGILITY OF QUANTUM CORRELATION TOWARD IMPERFECT DETECTION 

This section aims at identifying regimes suitable for high sensitivity image detection, which requires quantum 
correlation and a large number of photons per pixel at the same time. The latter condition is necessary because of 
the unavoidable presence of detection noise. If we denote by {SN'^)i,kg = {SN-^)bkg + {5N2)bkg tbe variance of the 
difference of background noise from symmetrical pixels of the CCD, we need to achieve {5N^)iikg/ {Ni + N2) ^ 1. As 
it was clearly illustrated in [3, [l^ , the task to achieve strong local sub-shot noise correlation with a large number 
of detected photons is experimentally demanding. We wish to investigate under which conditions the SNR obtained 
with the PDC source is able to beat the SQL defined in Eq. PH)) . In particular we shall emphasize the negative role 
of the excess noise present in the signal and idler beams, showing that the quantum nature of the correlation becomes 
fragile against unavoidable imperfections in the detection procedure. Because of this feature, the conditions of high 
parametric gain and high excess noise in the experiment described in would not suit our purpose. Indeed, the 
experimental data displayed a transition from quantum (sub-shot noise) to classical (above shot noise) correlation as 
the number of PDC photons per pixel was raised above 15-20, a regime where the intrinsic noise of the CCD is not 
negligible. 

It is well known that the signal and idler fields taken separately display a thermal-like statistics. The variance of 
the photon number intercepted by the pixel detectors can therefore be written in the form [2^ 



\2 



(SNi) = {N2) + ^ (20) 

where M denotes the degeneracy factor representing the number of spatial and temporal modes collected by the 
detectors. (Accurate measurements of the photon statistics in the high gain regime of PDC are described e.g. in f23|.) 
According to definition ([8]), it follows therefore that the excess noise coincides with the mean number of photons per 
mode: 

En^^. (21) 

This quantity depends only on the parametric gain and is on the order of sinh^ g. 

A large amount of excess noise has a detrimental effect on quantum correlation as one takes into account small 
unbalances between the test and the reference arms of the imaging system, which in real experiments can never be 
completely suppressed. Let assume for example that the signal and the idler fields are not perfectly balanced, but 
rather undergo different losses in the test and the reference arms, either during propagation, or in the detection 
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process. Indicating with ryi and 772 the different effective quantum efficiencies in the two detection channels, it is found 
that 



cr = 1 - 77 H — I E„ 

277 



(22) 



where fj = (771 +772)72. Thus we see from relation (j22p that a necessary condition to have sub-shot-noise correlation 
is that En ^ 277/(771 — 772)^. Because (771 — 772)^ is usually small compared to unity, the condition can be generally 
fulfilled, e.g. by applying a careful compensation of the losses in the two detection channels [lol [T3|. 

However, when dealing with an imaging scheme as in our case, another source of unbalance which is usually more 
difficult to control experimentally derives from the inaccuracy in the determination of the center of symmetry of the 
far-field pattern in the detection plane. Typically, the distance x shift between the selected and the exact centre of 
symmetry is on the order of one fourth the size of a CCD pixel (20/im in the experiment 14]) and is small compared 
with the far field coherence length Xcoh (a few tens of microns in [14]) but not negligible. In the following, all the 
quantities of interest will refer to a pair of pixel detectors located at symmetrical positions with respect to the optical 
axis as the one shown in Fig. [2l In particular A^i and N2 will indicate the operators associated with the number of 
photons measured by such a pair of detectors. 
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FIG. 4: Plot of the correlation factor ct as a function of x shift for increasing values of the pump pulse duration, considering 
a pixel size of (a) and 40/im (b). The value of the parametric gain (see legend) changes with Tp so that the number 

of detected photons is fixed, about 3500 ph. /pixel in case (a). The pump beam waist is Wp — 1500/.im, the lens focal length 
/ = 5cm, the detection quantum efficiency 77 = 0.9. 

Using the numerical model illustrated in Sec lIV| we investigated the behaviour of the correlation factor a toward 
errors in the center of symmetry Xshift, considering different values of the parametric gain g in order to vary the 
amount of excess noise. The far field measurement scheme considered in the simulation is identical to the one 
illustrated in Fig. ^ except we do not include the absorbing object, as we are presently interested only in evaluating 
the correlation factor. The results are illustrated in Fig. 21 We considered decreasing values of the parametric gain 
g in order to reduce progressively the amount of excess noise (which behaves as sinh^ g) starting from rather high 
values {En ^ 100 for g — 4.15) down to a small amount (i?„ ~ 0.5 for g = 1.45). At the same time the pump pulse 
duration Tp is increased as indicated in the figure legend so that the number of photons (A^i) detected by the pixel 
detectors remains unchanged, about 3500 ph. /pixel. In such conditions, the excess noise decreases as the inverse of 
the pump pulse duration Tpump- The two plots corresponds to two different sizes of the pixel detectors: 20/im in case 
(a) and 40/j,m in case (b). Referring to the 20/j,m pixels of the CCD used in experiment [13], the 40/Ltm pixel size can 
be obtained by performing a 2x2 pixel binning. With our choice of the lens focal length (/ = 5c77i) and the pump 
beam waist {wp — 1550/im), the typical coherence length in the detection plane is of the same order of magnitude of 
the pixel size or slightly smaller. Except for the case Tp = 5ps, the long duration of the pump pulse with respect to 
the PDC coherence time {Tcoh ~ Ips) did not allow to simulate PDC propagation using a unique numerical array. For 
this reason we divided the pump temporal profile in a sequence of 5 ps long intervals over which the generated PDC 
field can be considered as nearly stationary. We performed simulations for each of these time intervals separately, 
summing up their independent contributions in order to evaluate the overall photon statistics. 
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As can be inferred by comparing Fig. lUJa) and (b), under the same gain/pulse duration conditions the 2x2 pixel 
binning [Fig.lDJb)] allows to enhance the correlation substantially (as long as Xghift ^ Xcoh) since the detection area 
covers a larger portion of the coherence area. Another relevant issue which can be inferred from the both plots is the 
deterioration of the correlation that occurs at high parametric gains even for x shift = 0. This effect is related to the 
broadening of the signal-idler correlation function which occurs at high parametric gains as discussed in SecfVl 
It becomes particularly relevant for the 5ps pump pulse with g = 4.15 and the 20/im pixel detectors (black squares 
in Fig. IHa): the correlation factor a goes rather close to unity, i.e. close the shot noise level, since for such a high 
value of g the FWHM of the ri2, xfwhm ~ 23/im, becomes larger than the size of the pixel detectors. On the other 
side, considering the 40/im pixel detectors, we see that a approaches the optimal value 1 — 77 = 0.1 as the pump pulse 
duration is increased and the parametric gain is decreased. 

From the same plots we see that a starts from its minimum value taken at x shift — and increases with x shift 
until it saturates to the value a sat = 1 + En as x shift becomes larger than Xcoh- It can be easily verified that the 
saturation value asat corresponds to the limit in which the photon number fluctuations revealed on the two detectors 
becomes completely uncorrelated, as the error in the symmetry center becomes larger than the width of signal-idler 
correlation function (which can be identified with Xcoh)- Considering the worst case Tp = 5ps, the large parametric 
gain makes the excess noise extremely high (£'„ ~ 100). As a consequence, we see that the slightest inaccuracy in 
the position of the symmetry center leads to a transition of the degree of correlation from sub-shot noise to strongly 
above shot noise (black square in Fig. [¥]). As the parametric gain is decreased a saturates slightly faster (since Xcoh 
diminishes) but to lower and lower values which tend to approach the shot noise level as i?„ goes to zero (see stars 
of the lOOOps pump pulse case). As M cx Tp becomes larger than (A^i) the photon statistics on the two detectors 
becomes indeed nearly Poissonian, so that {5Nl) « (A^i) (see Eq. and a is at worse on the order of unity, the 

value corresponding to the standard quantum limit. 

If we approximate the curves of Fig. 2] with straight lines that go from the minimum value ^1 — 7] (at Xshift = 0) 
to the saturation value Gsat = 1 + En (at Xshift — Xcoh), we can obtain a rough estimate of the error shift value for 
which the transition from sub-shot-noise to above shot-noise correlations occurs, i.e. 



a > 1 



shift 
•^coh 



> 



1 + En/V 



(23) 



A more precise analytical derivation of the dependence of a on Xshi ft will be given in [24[ ; however we can infer from 
the behavior illustrated in Fig. 3] and from expression (|23p that the precision in the determination of the symmetry 
center becomes an extremely relevant issue from the experimental point of view whenever the excess noise is a large 
quantity, i.e. for short pump pulses and high parametric gains. Moreover, relation (j23p does not takes the broadening 
of the correlation function occurring with increasing g and therefore the fragility of the quantum correlation is 
underestimated. Clearly, a way to compensate for this negative effect would be to increase the detection area. On 
the contrary, when the excess noise is lowered by decreasing the parametric gain, the quantum correlation becomes 
more robust. 




x„=250ps 

SQL 
x-=1000ps 



6000 



FIG. 5: Plot of the correlation factor cr as a function of (A'^i) for different values of the pump pulse duration Tp (see legend). 
The error in the symmetry center is x shift = 4/im. (A''i) is increased by varying the parametric gain g. The pixel detector size 
is 20/xm. The other parameters are the same as in Fig. [l] 



The previous results clearly suggests that sub-shot noise correlation can be maintained at higher photon number 
values if the parametric gain is reduced and the pump pulse duration is increased, typically in the nanosecond range. 
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The simulation shown in Fig. [5] confirms this hypothesis: it plots a as function of the mean photon number per pixel 
(TVi) considering a fixed error in the determination of the symmetry center, x shift = 4 ^m, and for increasing values 
of the pump pulse duration, starting from 5ps up to Ins. Each point of the curves is obtained for different values 
of the parametric gain g in the the range between 1.2 and 3.5, while the solid lines correspond to a linear fit of the 
numerical data. As a general behaviour, we see from the figure that a increases linearly with (A^i) with a slope which 
decreases as the pump pulse duration is increased. We verified that the angular coefficient of the linear fits (solid 
lines) scales as the inverse of M ~ Tp/Tcoh- For the Ins pulse one can therefore reach sub-shot noise correlation at 
much higher photon number values than for shorter pulses. 



VII. SIMULATION OF THE HIGH-SENSITIVITY IMAGING EXPERIMENT 



From the previous results we have seen that in the presence of excess noise a careful balance of the two detection 
channels and a precise determination of the symmetry center are necessary in order to approach the (ideal) minimum 
value a — 1 ~ rj. On the other hand, we see both from Eq. and relation that this requirement becomes less 
stringent as the excess noise is reduced by lowering the parametric gain, so that M becomes on the same order of 
magnitude of (A^i) or larger. We can take advantage of this feature by considering long pump pulses, typically in the 
nanoseconds range, so that Tp ^ T^oh Etnd a large number of temporal modes are amplified at a low parametric gain, 
keeping large the total number of photons collected by the pixel detectors. 
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FIG. 6: Evaluation of the SNR obtained by increasing the pump pulse duration for a fixed parametric gain g — 1.45. The 
object absorption coefficient is a = 0.04. In (a) the pixel size is 20/im, in (b) a 2x2 binning is performed giving a pixel size 
of 40/im. The lower horizontal scale gives the number of detected photons per pixel, the upper scale the corresponding pump 
pulse duration. The evaluated degree of correlation a and the SNR improvement R axe indicated. The test and reference arms 
are assumed to be perfectly balanced (in particular x shift = 0). The other parameters are the same as in Fig. H) 




FIG. 7: Simulation of the retrieval of a weak object (a) through a diflferential measurement with a 125ps pump pulse and 
~ 1750 ph. /pixel. The indicated SNR values corresponds to those with the PDC source (b) and with the classical source (c) in 
the simulation of Fig. |6jb) with the 40/im pixel size (see dotted vertical line in this figure) . 

We now illustrate the numerical simulation of a high-sensitivity experiment performed considering the measurement 
scheme of Fig. O Let us first consider the case in which test and reference detection channels are perfectly balanced 
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(in particular x shift = 0). The object we considered took the form of a weakly absorbing mask on the numerical grid 
with a constant absorption coefScient a — 0.04, as the one shown in Fig.[7Ka). FigurelHlillustrates the behaviour of the 
SNR as a function of pump pulse duration. In those simulations the parametric gain is kept constant [g = 1.45), while 
Tp is increased up to lOOOps with the mibinned pixels [Fig. [Sfa)], up to 250ps with the 2x2 pixel binning [Fig. [BUb)]; 
in such a way (A'^i) (lower horizontal scale) increases almost linearly with the pump pulse duration (upper axis) up 
to 3500 photons per pixel and the SNR increases with the square root of (A^i) as predicted by Eqs. ([9]) and ((7]). 

The lower curves (black squares) represent the SQL obtained with the splitted coherent beams, while the upper 
curves (hollow triangles) reproduce the results obtained with the PDC source. With the chosen parameters, the 
coherence length ojfwhm in the detection plane was evaluated on the order of 10/xm. As discussed in SeclVIl (see 
Fig. m, only the 40/xm pixels size allows to approach the optimal degree of correlation determined by the finite 
quantum efficiency of the detectors, i.e. cr = 1 — 77 = 0.1. For this reason the factor of improvement R obtained with 
the unbinned pixels [Fig. Elja)] is only ^ 1.9, while with the 2x2 pixel binning [Fig. HKb)] we are close to the optimal 
value R = 2.9 (obtained from Eq. pT|) with a = 0.1, En, « 0.5). From these results it clearly emerges that the size of 
the coherence area determined by the optical setup imposes a lower limit to the size of the object details that can be 
resolved with a sensitivity beyond the SQL. 

The dashed lines in the figure are obtained by using the expressions of the two-mode calculations ([7]) and (O and 
fit well the numerical results. The values of the relevant fitting parameter a are indicated beside each plots (£"„ is on 
the order of unity and does not affect the SNR value significantly, as can be inferred from expression ([7])). It should 
be noticed that both R and a are almost independent on the pump pulse duration, as long as the parametric gain g 
is fixed. 
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FIG. 8: The ratio R — SNRct/SNRsql is plotted as a function of x shift for different values of the pump pulse duration, 
considering 20/im (a) and (b) pixel detectors. All other parameters are the same as in Fig. 3] 

The images of Fig. [7] allow to appreciate visually the improvement that can be achieved with the PDC source (b) 
with respect to the classical source (c) with the 40/im pixel size. The indicated SNR values are those obtained with 
a 125ps pump pulse (giving ^ 1750 photons/pixel) corresponding to the dotted vertical line in Fig. [HI 

Finally we investigated on the sensitivity of the SNR improvement with respect to inaccuracies in the determination 
of the center of symmetry of the PDC far field. In Fig. [8] we plot the ratio R = SNR^/ SNRsql as a function of 
X shift for different values of the pump pulse duration and for a pixel detector size of 20^m (a) and 40/Ltm (b). The 
chosen parameters are the same as in Fig. |4j the parametric gain g is diminished for increasing values of Tp, so that 
the number of photons detected on the pixel area remains unchanged, about 3500 ph. /pixel. 

A first relevant issue is the deterioration of the correlation and the SNR that occurs at high parametric gains even 
for Xshift = 0. This effect is related to the broadening of the coherence area which occurs at high parametric gains 
due to the finite transverse size of the pump pulse discussed in|Vl It becomes particularly relevant for the 5ps pump 
pulse with g ~ 4.15 (black squares in the figures): the correlation factor a goes rather close to unity (see Fig. 2]) and 
the SNR is at best close to the SQL value, for Xshift — 0, as can be seen from Fig. [8l 

Furthermore, for Tp = 5ps the excess noise is extremely high {En ~ 100). As a consequence, we see that the slightest 
inaccuracy in the position of the symmetry center leads to a transition from sub-shot noise to strongly above shot 
noise signal/idler correlation (black square in Fig. [5]). For this reason the SNR rapidly decreases well below the SQL 
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as X shift increases, as can be seen from Fig. |8l The situation improves as the pulse duration is increased and the 
parametric gain is decreased, since the excess noise diminishes and the SNR becomes more robust toward Xsm ft ■ For 
example, considering the best detection conditions with the 40/Ltm pixel size [Fig. [8)3] , the tolerance on XsM ft (in order 
to have i? > 1) is less than ^ 2/im for g — 4.15, while it increases to ^ 14/im with g = 1.45. 

To conclude, our simulations show that the detrimental effects of the broadening of the coherence area added to 
the fragility against imperfections in the measurement process become a much less relevant issue in the regime of 
low-gain, long-pulse regime. 



We have proposed an imaging scheme that uses the multi-mode twin beams generated in the process of single-pass 
parametric down-conversion to detect the spatial distribution of faint objects through a differential measurement 
technique. The object is located in the PDC far-field and the image is retrieved by subtracting pixel by pixel the 
measured signal and idler photon fluxes. In this way the detection scheme can achieve a sensitivity beyond the SQL 
by exploiting the local character of the quantum correlations of the PDC far field, a property which was demonstrated 
experimentally in (author?) [l3. [l5| in a regime of high parametric gain. 

We developed a stochastic numerical model which includes realistic features of the imaging system in order to verify 
to which extent the PDC source allows to enhance the SNR of the measurement with respect to that obtained from a 
classical source. We pointed out the fragility of the sub-shot noise correlation against imperfections of the detection 
system in presence of a large amount of excess noise. In particular we analyzed the effect of small errors (on the order 
of a few microns) in the determination of the center of symmetry of the PDC far field, showing that they lead to a loss 
of correlation proportional to the excess noise of the source. Other experimental imperfections, such as unbalances 
due to different losses in the test and reference arms and light scattering from the background, produce a similar 
deterioration of the sub-shot noise correlation. The quantum character of the correlation is therefore rapidly lost as 
the parametric gain (and therefore the excess noise) is increased. 

We suggest that the use of pump pulses with duration on the order of a few nanoseconds or more, much longer than 
the PDC coherence time, would allow to generate the same number of photons than in the case of the picoseconds 
pulses used in [l4| with a much lower parametric gain and a negligible excess noise. Our simulations demonstrate that 
this represents a clear advantage in that the sub-shot noise correlation becomes more robust and can be maintained 
up to gains corresponding to several thousands of PDC photons per pixel, a regime where the noise of the CCD can 
be neglected [25j. We used the same numerical model to simulate the imaging experiment, showing that under these 
conditions the use of a PDC source lead to a substantial SNR enhancement with respect to classical imaging in a 
domain of parameters accessible to experimental implementations. 

Moreover, from our analysis it emerges that the degree of quantum correlations, and hence the improvement in the 
SNR depends critically on the ratio between the pixel area and the far-field coherence area. The latter sets a lower 
limit on the size of the details that can be resolved with sensitivity above the SQL, and thus provides a resolution 
limit for high-sensitivity imaging with twin beams. 



This work was carried out in the framework of the PRIN project of MIUR "Twin beams in Quantum Imaging 
applications and metrology" . 
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